Deformation effect on total reaction cross sections for neutron-rich Ne-isotopes 
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Isotope-dependence of measured reaction cross sections in scattering of 28_32 Ne isotopes from 12 C target 
at 240 MeV/nucleon is analyzed by the double-folding model with the Melbourne g-matrix. The density of 
projectile is calculated by the mean-field model with the deformed Wood-Saxon potential. The deformation 
is evaluated by the antisymmetrized molecular dynamics. The deformation of projectile enhances calculated 
reaction cross sections to the measured values. 
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I. INTRODUCTION 

Exploring unstable nuclei is one of the most important sub- 
jects in nuclear physics. Actually, it was reported that unstable 
nuclei have exotic properties such as the halo structure lHJ-Ht] 
and the loss of magicity for nuclei in the so-called "Island of 
inversion". The term "Island of inversion" was first introduced 
by Warburton Q4[| to the region of unstable nuclei from 30 Ne 
to 34 Mg. In the region, the low excitation energies and the 
large B(E2) values of the first excited states suggest strong 
deformations @-03], which indicates that the N = 20 magic 
number is no longer valid. These novel quantum properties 
have inspired extensive experimental and theoretical studies. 

Important experimental tools for exploring unstable nuclei 
are the reaction cross section ctr or the interaction cross sec- 
tion o~i and the nucleon-removal cross section ct_„ with ra- 
dioactive beams llT-S [loll : for the scattering of unstable nu- 
clei, o~i agrees with or in general, since projectile excitations 
to its discrete excited states do not exist. Very recently, ctt 
was measured by Takechi et al. ifTHl for 28_32 Ne located near 
or in "Island of inversion". Furthermore, a halo structure of 
31 Ne was reported by the experiment on the one-neutron re- 
moval reaction fl3l . This is the heaviest halo nucleus in the 
present stage suggested experimentally and also reside within 
the region of "Island of inversion". 

As a useful theoretical tool of analyzing or, we can con- 
sider the microscopic optical potential constructed by the 
double-folding model (DFM) with the q-matrix effective 
nucleon-nucleon (NN) interaction Ill3l42lll . when the pro- 
jectile breakup is weak. For the nucleon-nucleus scatter- 
ing, the single-folding model with the g-matrix well repro- 
duce the data on or and the elastic-scattering cross section 
systematically lfl3ll . For the 31 Ne scattering from 12 C at 
240 MeV/nucleon, the breakup cross section is at most 1.5% 
of or l22tl . Hence, DFM is applicable also for analyses 
of measured isotope-dependence of or in the scattering of 
28-32 Ne from 12 C target at 240 MeV/nucleon ifTTfl . 

In DFM, the g-matrix is folded with the projectile and tar- 
get densities. If the projectile deforms, the density profile 
changes; the surface diffuseness increases because of the elon- 
gation. This gives rise to the effective growth of the root- 
mean-square (RMS) radius and eventually the increase of or. 
Therefore, the amount of deformation is important. Nuclei in 



the island of inversion are spherical or only weakly deformed 
in the Sky rme and/or Gogny HF (HFB) calculations; see, e.g., 
Refs. J23ll24ll . It is even pointed out that the observed large 
B(E2; 2 + — > + ) values can be understood as a large ampli- 
tude vibration around the spherical shape rf25ll . In such a sit- 
uation, the additional correlations by the angular momentum 
projection (AMP) often leads to possible deformed shapes; 
see Ref. 112611 for Ne isotopes. 

Recently a systematic investigation employing the antisym- 
metrized molecular dynamics (AMD) with the Gogny D 1 S in- 
teraction has been performed for both even and odd N nuclei 
in the island of inversion Ir27ll . The AMD (with AMP per- 
formed) gives rather large deformations, which is consistent 
with the AMP-HFB calculations (UHal- A consistent pic- 
ture of even and odd isotopes has been obtained, where the 
n-particle m-hole excitations of the Nilsson orbits play im- 
portant roles to determine deformed configurations. Although 
it is difficult to distinguish the dynamic shape-fluctuation and 
static deformation in these light mass nuclei, one may use the 
deformed shape suggested by the AMD calculation to see its 
effect on or. 

In this paper, we analyze the measured isotope-dependence 
of ctr in scattering of 28 ~ 32 Ne isotopes from 12 C target 
at 240 MeV/nucleon, using DFM with the Melbourne g- 
matrix f 13] and the deformed projectile density suggested by 
the AMD calculation. 



II. THEORETICAL FRAMEWORK 

A microscopic optical potential U between a projectile (P) 
and a target (T) is constructed with DFM. The direct and ex- 
change parts, Uu and Uex, are obtained by ll28ll29ll 

Ud(R) = J pp(r P )p T (r T )v D (p, s)dr P dr T , (1) 

Uex(R) = J Pp(rp,r P + s)p T (r T ,r T - s) 

x vex(p, s) cxp [iK(R) ■ s/M]drpdr T , (2) 

where s = rp — + R for a position vector R of P from 
T. The original form of Uex is a non-local function of R, but 
it has been localized in Eq. (O with the local semi-classical 
approximation (3, where HK(R) is the local momentum 



of the scattering considered and M — ApA-r/ (Ap + At) 
for the mass number Ap (At) of P (T). The validity of this 
localization is shown in Ref. l30ll . Here, the effective NN 
interactions, vp> and «ex> ar e assumed to depend on the local 
density 



p = pp(rp + s/2) + pt (t*t - s/2) 



(3) 



at the midpoint of the interacting nucleon pair. 

The microscopic potential U is not spherical, if one or both 
of the densities pp and pt are non-spherical. As shown in 
Ref. El, however, the effect is found to be negligible for 
heavy-ion scattering. For intermediate incident energies of 
our interest, this can be understood with reasonable approxi- 
mations. For such energies, the rotational motion of deformed 
P (rotor) is negligible compared with the center-of-mass mo- 
tion of P. Hence, the adiabatic approximation is applicable for 
the rotational motion. Under the approximation, P is scattered 
by U depending on not only R but also SI the direction (Euler 
angles) of the symmetry axis of rotor in the space-fixed frame. 
Using the eikonal approximation for the center-of-mass mo- 
tion of P, one can obtain 



°R = J db(l - \S\ 2 ) 



with 



S 



dn 

8^ 



I 

hv 



dzU(R, n) 



(4) 



(5) 



for R = (b, z) and the velocity v of P. Here, we have as- 
sumed that the ground state of P is a + state for simplicity. 
The S'-matrix element S includes back-coupling effects of the 
rotational excitations on the elastic scattering. Separating U 
into the spherical and non-spherical parts, Uq and AU, we can 
get 



with 



S — Sn + Sn 
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(7) 

(8) 



The non-spherical correction to the spherical part So starts 
with 5 2 , since the correction of order S vanishes because of 
the angle average. The leading-order correction is significant 
only at large b, since Sq vanishes at small & as a result of 
the strong absorption. For large b, 5 keeps small because v 
is large and the range of the z integration is small. Actually, 
we confirmed through numerical calculations that the leading- 
order correction to or is 0.01% for the 30 Ne+ 12 C scattering at 
240 MeV/nucleon; note that the error of the eikonal approx- 
imation is less than 1% for this scattering. Also for the case 
that the spin of P in its ground state is non-zero, it is possible 
to prove that the leading-order correction to or is of order S 2 . 



Thus, the effect of AU on or is negligible. Therefore, only 
the spherical part of the density is taken in this paper. De- 
tailed discussion on this non-spherical effect will be made in 
the forthcoming paper. 

As for pt, we use the phenomenological 12 C-density de- 
duced from the electron scattering ll3~i1l by unfolding the finite- 
size effect of the proton charge in the standard manner 1I32I . 
Meanwhile, pp is calculated by the mean-field model with 
a given average potential or with the self-consistently deter- 
mined potential by the Hartree-Fock (HF) method. No effect 
of pairing is included for simplicity. The Ne isotopes (pro- 
jectiles) under discussions are supposed to be in the island of 
inversion (or at its boundary), and expected to be strongly de- 
formed. In order to investigate the effect of deformation, we 
take a deformed Woods-Saxon (WS) potential i33ll . in which 
the axially deformed surface S(/3) is specified by the radius, 



R(0;/3) = R c v (p)[l + E\P\Yxo(9)]> 



(9) 



with the deformation parameters (3 = {(3\} and a volume 
conserving factor c v (/3). The potential value is determined by 
replacing the quantity (r — Rq) in a spherical potential to the 
distance from the surface S(/3) (with minus sign if the point is 
inside it). The Coulomb potential created by charge (Z — l)e 
distributed uniformly inside the surface S(/3) in Eq. (O is in- 
cluded for protons. The single-particle eigenstates are cal- 
culated by the (cylindrical) harmonic oscillator basis expan- 
sion. More than twenty oscillator shells are included and the 
convergence of the result is carefully checked to obtain reli- 
able density distributions. The nucleon density is obtained by 
summing up the contributions of occupied Nilsson levels. The 
density distribution thus calculated pp (r, 6) is the one in the 
intrinsic (body-fixed) frame, and depends on the polar angle 
9 from the symmetry axis. As mentioned above, the density 
in the laboratory frame used in DFM (Eqs. ([TJ and (f2}) is ob- 
tained by the angle average: 



1 f 77 

Pp(r) = 2 y o P P n) (r,6»)sin 



(10) 



We have checked that the angle-averaged density agrees with 
high accuracy with the density calculated by the angular mo- 
mentum projection from the Slater determinantial wave func- 
tion composed of the occupied WS orbits. 

No center of mass (CM) correction is included for the cal- 
culation of the density. We have checked by the spherical 
Gogny HF calculation that the CM correction (including the 
two-body contributions) to the RMS radius is about 1% reduc- 
tion for all the isotopes. The amount of reduction is smaller 
than the enhancement caused by the deformation effect, but is 
non negligible; we will return to this point latter. 



III. RESULTS 

We test the accuracy of DFM with the Melbourne g-matrix 
for 12 C+ 12 C scattering at 250.8 MeV/nucleon. As shown in 
Table H or calculated with the Melbourne g-matrix is con- 
sistent with the experimental data; more precisely, the latter 



TABLE I: Reaction cross sections for C+ C scattering at 
250.8 MeV/nucleon for two types of effective nucleon-nucleon in- 
teractions. The cross sections are presented in units of mb. 



Exp. [34] Love-Franey r351 Melbourne- g f 131 
782.0 ± 10 918 796 



TABLE II: Deformation parameter P2 used in the calculation of 
density of Ne isotope (those with higher multipoles A > 2 are not 
included). 



nuclide 


28 Ne 29 Ne 


30 Ne 31 Ne 


32 Ne 


ft 


-0.291 0.445 


0.400 0.422 


0.335 



is slightly smaller than the former by the factor F = 0.982. 
The table also shows the result of the Love-Franey /-matrix 
nucleon-nucleon interaction in which the nuclear medium ef- 
fect is not included. The difference between the two theo- 
retical results is about 122 mb that corresponds to 16% of 
the experimental data. Thus, the medium effect is impor- 
tant at this incident energy. For the 27 A1 + 12 C scattering 
at 250.7 MeV/nucleon, or calculated with the phenomeno- 
logical 27 Al density [31] and the normalization factor F is 
1 164 mb, while the experimental value is 1159 ± 14 mb l34tl . 
The normalization procedure thus justified is applied for the 
28-32j s j e + i2£ sca ttering at 240 MeV/nucleon analyzed be- 
low. 
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FIG. 1: (Color online) Reaction cross sections for scattering of Ne 
isotopes from 12 C at 240 MeV/nucleon. The results obtained with 
the calculated density with the Woods-Saxon potential are denoted 
by dashed line, and those with the Gogny HF by solid line. The 
spherical shape is imposed. The nuclei with A > 30 are unbound. 
The experimental data are taken from Ref. iTTTll . 



As for the parameter set of the WS potential, i.e., the depth, 
radius and the diffuseness of the central as well as the spin- 
orbit potentials, we employ the one provided recently by 
R. Wyss HH; see Table I of Ref. l37ll for the actual values 
of parameters. This set is intended to reproduce the spectro- 
scopic properties of high-spin states from light to heavy de- 
formed nuclei, e.g., the quadrupole moments and the moments 
of inertia, and at the same time the RMS radii crucial for the 



present analysis. In order to check that the present WS poten- 
tial gives reasonable results, we compare in Fig.[T]the reaction 
cross sections calculated by using two densities; one obtained 
by the Gogny D1S HF calculation and another with the WS 
potential: The spherical shape is imposed with the filling ap- 
proximation in this calculation. A good agreement shown in 
the figure indicates that the density distributions in the two 
models are similar, which is also confirmed by the calculated 
RMS radii (see Fig. EJ. 

The reaction cross section is sensitive to the amount of de- 
formation. We then employ the deformed shapes suggested 
by the AMD calculation to see the effect on or. As a sim- 
ple estimate we only include the Y20 deformation in Eq. (0 
and the deformation parameter $2 in each isotope is deter- 
mined to reproduce the calculated ratio of RMS radii along 
the long and short axes by AMD; the resultant values used in 
the following analyses are given in Table [XT] With these P2 
values, the Nilsson orbits of the last-odd-neutron in 29 Ne and 
31 Ne are [200]l/2 and [321]3/2, respectively, in accordance 
with the AMD calculation. Note that the nucleus 28 Ne is at 
the boundary of the "Island of inversion", and AMD predicts 
strong mixing between the states with oblate and prolate de- 
formation. In the present calculation, we have employed the 
/?2 value of the oblate minimum, which is the main compo- 
nent. 
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FIG. 2: (Color online) Reaction cross sections for scattering of Ne 
isotopes from C at 240 MeV/nucleon. The dashed and solid lines 
represent results of the spherical and deformed WS pote ntials, re- 
spectively. The experimental data are taken from Ref. in. 



The results of or including the effect of quadrupole defor- 



mation (see Table HB are shown in Fig. [2] Compared to the 
results with the density of the spherical cases, the effect of 
deformation increases the cross section considerably. The en- 
hancement makes the calculated cross sections almost consis- 
tent with the measured cross sections for 28_32 Ne, although 
the difference of or between 30 Ne and 31 Ne is small in the 
model calculation compared with the difference deduced from 
the measured cross sections. We will return to this point later. 
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to that of or in Fig. [2 which is reasonable because of a simple 
estimate, 

a R « 7 r[(r 2 )J ) /2 + <r 2 )^ 2 ] 2 , (11) 

where (r 2 )^/ 2 is the RMS radius for target. Note that the 
amount of increase of the RMS radii from the spherical shape, 
which is roughly proportional to /3f > is only 4 — 6%: It is sur- 
prising that such a small effect is detectable in experimental 
data. The present analysis clearly tells us that or reflects very 
precise information, and its measurement is extremely useful 
to study the nuclear structure of unstable nuclei. The radial 
dependence of the matter density is plotted in Fig. |4] The 
deformed WS density (solid curve) is enhanced by the defor- 
mation effect from the spherical WS density (dashed curve) 
at r > 4 fm. The enhancement of the reaction cross section 
is caused by that of the density in this tail region. This is the 
main reason why we do not directly use the calculated AMD 
density, which decreases more rapidly in the tail region be- 
cause of the usage of one-range gaussian wave functions. 
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FIG. 3: (Color online) Matter RMS radii for Ne isotopes for 
the spherical WS potential (crosses), the deformed WS potential 
(squares) and the Gogny HF (circles). 
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FIG. 4: (Color online) The radial dependence of the matter density 
for 29 Ne. The dashed and solid lines show results of the spherical 
and deformed WS potentials, respectively. 

The increase of or caused by the deformation can be rather 
nicely understood if one looks into the (matter) RMS radii 
(r 2 )p shown in Fig. [3] They are calculated by using the pro- 
jectile density pp based on the spherical and deformed WS po- 
tentials. The increase of (r 2 )p 2 in Fig.[3]nicely corresponds 



IV. DISCUSSIONS 

The enhancement of the reaction cross sections caused by 
the deformation effect is conspicuous as shown in Fig. [2] The 
enhancement makes the calculated cross sections almost con- 
sistent with the observed ones for 28 ~ 32 Ne; more precisely, 
the calculated cross sections slightly overshoot the data for 
30 ' 32 Ne, but slightly undershoot the data for 31 Ne. In the 
spherical HF calculation, the CM correction to the RMS radii 
yields 1% reduction. This leads to 1.1% reduction of ctr 
through relation (fTTT i. It is very likely that the deformed WS 
model has almost a same amount of or reduction. After this 
reduction, the calculated cross sections agrees with the data 
for 2 8-30,32 Ne ^ but underestimates the data by 32 ± 22 mb for 
31 Ne. Thus, the theoretical results are consistent with the data 
for 30 Ne but not for 31 Ne, because the difference of or be- 
tween the two nuclei is smaller in the model calculation than 
in the data. 

The difference of or between 31 Ne and 30 Ne corresponds 
to the one-neutron removal cross section of 31 Ne, if the 
breakup cross section of 31 Ne is negligible l22"tl . The differ- 
ence between the observed reaction cross sections is 86 mb, 
while the direct measurement on the one-neutron removal 
cross section yields 79 mb fl2Tl . Thus, the two experimental 
data are consistent with each other, indicating that the breakup 
cross section is small. Meanwhile, the difference of the calcu- 
lated reaction cross sections between the two nuclei is 41 mb 
and smaller than the experimental results. 

As for 31 Ne, the single-particle energies of the last neutron 
are about —2 MeV in the present deformed WS potential with 
fii value given in Table HI1 The underestimation of the present 
value for 31 Ne may mean that either the depth of the present 
WS potential is too deep or @2 is too small. For example, 
compared with the WS potential in Ref. ll38tl . the binding en- 
ergies of relevant Nilsson orbits are about 2 MeV larger in the 
present case, though the Nilsson diagrams are very similar to 



each other. It turns out that we can obtain good agreements of 
or for 31 Ne either by shallowing the potential depth by factor 
0.943 or by increasing the deformation up to fa = 0.590. 

In the case of 31 Ne, its spin-parity and neutron config- 
uration are still under debate. Our prediction of the last- 
odd-neutron orbit is [321]3/2 with the single-particle energy 
— 1.947 MeV. The energy increases to —0.974 MeV when the 
potential is reduced by a factor 0.943 to account for the ob- 
served central value of or, while the last-odd-neutron orbit 
changes to [200] 1/2 and the energy decreases to —2.803 MeV 
when fa is increased to 0.590. The measured separation en- 
ergy of 31 Ne, 0.29 ± 1.64 MeV (H], is more consistent with 
the single-particle energy of the shallower potential rather 
than that of larger fa. 

It should be mentioned that the present calculation of ctr is 
not sensitive to the isovector properties, e.g., the neutron skin. 
Although the matter radii calculated with the present WS and 
with the Gogny D1S HF (imposing the spherical shape) per- 
fectly agree and so do the reaction cross sections (see Fig.Q~|i, 
the skin thicknesses in the two calculations are rather differ- 
ent: e.g., (r 2 ), 1 / 2 - (r 2 )p /2 w 0.67 and 0.41 fm with the WS 
and the Gogny HF, respectively, in 30 Ne. Additional informa- 
tion is necessary to probe the property like the skin thickness. 

V. SUMMARY 

Isotope-dependence of measured reaction cross sections 
in scattering of 28_32 Ne isotopes from 12 C target at 



240 MeV/nucleon is analyzed by the double-folding model 
with the Melbourne g-matrix. The density of projectile is cal- 
culated by the mean-field model with the deformed Wood- 
Saxon potential. The deformation is evaluated by the antisym- 
metrized molecular dynamics. The deformation of projectile 
enhances calculated reaction cross sections to the measured 
values. The increase of the RMS radii by the deformation is 
only 4 — 6%, but it is quite important that such a small effect 
is detectable in the experimental data. Owing to this effect, 
the calculated reaction cross sections reproduce the data for 
28-30 ' 32 Ne. For 31 Ne, however, the present results still un- 
derestimate the measured cross sections. The underestimation 
may suggest that the extra weak-binding effect for neutrons 
plays an important role particularly for 31 Ne. 
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